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ABSTRACT 

We construct models of core-collapse supernovae in one spatial dimension, including rotation, angular momen- 
tum transport, and viscous dissipation employing an a-prescription. We compare the evolution of a fiducial 1 1 Mq 
non-rotating progenitor with its evolution including a wide range of imposed initial rotation profiles ( 1 .25 < f o < 8 
s, where Pq is the initial, approximately solid-body, rotation period of the iron core). This range of Pq covers the 
region of parameter space from where rotation begins to modify the dynamics (Pq ~ 8 s) to where angular ve- 
locities at collapse approach Keplerian (Pq ^ 1 s). Assuming strict angular momentum conservation, all models 
in this range leave behind neutron stars with spin periods < 10 ms, shorter than those of most radio pulsars, but 
similar to those expected theoretically for magnetars at birth. 

A fraction of the gravitational binding energy of collapse is stored in the free energy of differential rotation. 
This energy source may be tapped by viscous processes, providing a mechanism for energy deposition that is not 
strongly coupled to the mass accretion rate through the stalled supernova shock. This effect yields qualitatively 
new dynamics in models of supernovae. We explore several potential mechanisms for viscosity in the core- 
collapse environment: neutrino viscosity, turbulent viscosity caused by the magnetorotational instability (MRI), 
and turbulent viscosity by entropy- and composition-gradient-driven convection. We argue that the MRI is the 
most effective. We find that for rotation periods in the range Po ^ 5 s, and a range of viscous stresses, that the 
post-bounce dynamics is significantly affected by the inclusion of this extra energy deposition mechanism; in 
several cases we obtain strong supernova explosions. 

Subject headings: stars: magnetic fields — stars: neutron — supernovae: general — hydromagnetics 



L INTRODUCTION 



Employing the now-standard suite of microphysics and a full 
solution to the Boltzmann equation for all neutrino species, the 
neutrino mechanism of core-collapse supernovae fails in one 
spatial dimension (Rampp & Janka 2000, 2002; Liebendorfer 
et al. 2001a,b; Mezzacappa et al. 2001; Thompson et al. 2003; 
Liebendorfer et al. 2004). The initial "bounce shock," formed 
as the core density reaches ^2 x 10'"* g cm"-* and nuclei undergo 
a phase transition to free nucleons, stalls almost immediately 
due to a combination of neutrino losses, the ram pressure of 
the infalling mantle, and the dissociation of nuclei into neu- 
trons and protons behind the shock. The revival of the shock to 
an asymptotic energy of order 10^' erg is the focus of modern 
supernova theory. Any mechanism for this revival must fun- 
damentally rely on the transfer of gravitational binding energy 
to the post-shock mantle. As its name makes clear, the "neu- 
trino mechanism" (Bethe & Wilson 1985) relies on neutrino 
interactions (primarily ^ pe~ and D^p ^ ne^) to transfer 
binding energy to the shocked matter in the post-bounce epoch. 
In spherical models, the stalled bounce shock remains trapped 
forever and continues to accrete the overlying progenitor Al- 
though some models show that multi-dimensional effects might 
be necessary for success of the neutrino mechanism (Herant 
et al. 1994; Burrows et al. 1995; Janka & Muller 1995, 1996; 
Fryer et al. 1999; Fryer & Warren 2002), multi-dimensional 
models employing more sophisticated neutrino transport fail 

' Hubble Fellow 

^ Astronomy Department & Theoretical Astrophysics Center, 601 Campbell Hall, 
^ Astronomy Department & Theoretical Astrophysics Center, 601 Campbell Hall, 
'* Steward Observatory, The University of Arizona, Tucson, AZ 85721; buiTows® 



to explode (Janka et al. 2002; Buras et al. 2003) - albeit 
marginally. 

In constructing one-dimensional models, the focus of the su- 
pernova community has been largely on refining details of the 
neutrino mechanism. These details include a fuller description 
of the nuclear physics and the neutrino interactions and, over 
the last ten years, a precise handling of neutrino transport (Mez- 
zacappa & Bruenn 1993a,b,c; Burrows et al. 2000; Rampp & 
Janka 2000, 2002; Liebendorfer et al. 2001a,b; Mezzacappa et 
al. 2001; Thompson et al. 2003). Precision neutrino transport 
is important because the degree of neutrino heating behind the 
shock controls the subsequent dynamics and because the opti- 
cal depth of this region is < 0.1. For this reason, the most re- 
cent one-dimensional models are fully spectral and multi-angle 
(Rampp & Janka 2000, 2002; Liebendorfer et al. 2001a,b; Mez- 
zacappa et al. 2001; Thompson et al. 2003); neutrino transport 
also presents the most significant technical hurdle in calculating 
realistic models of core collapse in two or three spatial dimen- 
sions (Livne et al. 2003). 

It is possible that yet more refined treatments of the radia- 
tion hydrodynamics or subtle changes in the microphysics of 
neutrino interactions may yield explosions in one-dimensional 
models via the neutrino mechanism. However, it is also pos- 
sible that the solution to the supernova problem does not lie 
in such technical improvements or in an accumulation of five 
and ten percent effects, but instead in new physics and order 
unity phenomena not previously identified or fleshed out. Such 
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improvements might include new massive pre-coUapse progen- 
itors with very different core entropy and lepton fraction, new 
finite temperature high-density nuclear equations of state (in- 
cluding exotic meson species, quark matter, phase transitions, 
etc. and the associated neutrino microphysics), or the interplay 
between rotation and magnetic fields. In this paper, we take a 
step towards addressing the last of these issues in detail. 

When the core of a rotating massive star collapses, some of 
the gravitational binding energy of collapse is stored in a strong 
shear profile. This free energy reservoir may be tapped by vis- 
cous processes. In this way, the binding energy of collapse, 
temporarily trapped in differential rotation, may be transferred 
to the matter as thermal energy. This is qualitatively different 
from the standard neutrino mechanism, which relies solely on 
neutrino interactions to communicate the binding energy of the 
collapsed iron core to the rest of the star We find that dissi- 
pation of shear energy can significantly affect the dynamics of 
supernovae; models that would otherwise fail succeed as a re- 
sult of viscous dissipation. The most promising source of vis- 
cosity is magnetic stresses generated by the magnetorotational 
instability (MRI; Akiyama et al. 2003) and perhaps magneto- 
convection (Duncan & Thompson 1992; Thompson & Duncan 
1993). Although magnetic fields in excess of 10^^ G are gener- 
ated, no large iron core seed field is required for the success of 
this mechanism, since any seed field grows exponentially on a 
rotation period. 

In §2, we provide some context and motivation for the 
present study and we review some previous research on core 
collapse with rotation. In §3, we describe the technical specifi- 
cations of our one-dimensional models including rotation, vis- 
cous dissipation of shear energy, and the corresponding trans- 
port of angular momentum. Section 4 describes mechanisms 
for viscosity, including neutrino interactions, the MRI, and hy- 
drodynamic convection. In §5, we present results from rotat- 
ing models of core collapse including the centrifugal force, but 
without viscous dissipation. We also provide a general discus- 
sion of the conditions for explosion and an explanation for why 
these models and non-rotating models fail. We compare the 
magnitude of the viscous processes discussed in §4 in §6. In §7, 
we present results from calculations of core collapse including 
viscous dissipation. We show that explosions are obtained, we 
explain why these models succeed, and we delineate the suc- 
cessful parameter space. In §8, we summarize our results and 
discuss the implications of our findings. 

2. SCENARIO & MOTIVATION 

When the iron core of a massive star collapses as a re- 
sult of the Chandrasekhar instabihty, the implosion is reversed 
at nuclear densities. Here, nuclei dissociate into free nucle- 
ons and the equation of state stiffens dramatically, driving the 
bounce shock wave into the supersonic infalling outer core. 
The bounce shock stalls almost iimnediately and a character- 
istic post-bounce quasi- steady- state accretion structure is ob- 
tained. The hot (temperatures of order 10 MeV) newly-born 
protoneutron star (PNS) has a neutrinosphere radius, Ri,, of ap- 
proximately 50 - 80 km (defined here by the radius at which 
the optical depth to electron neutrinos ^2/3, a very energy- 
dependent quantity). Overlying the PNS is a subsonic accretion 
flow, bounded by the stalled stand-off bounce shock at a radius 
of ~ 150 - 200 km. The region between the PNS and the shock 
consists of a cooling layer just above the PNS and a region of 



net heating (the "gain" region) between the cooling layer and 
the shock. The boundary between net cooling and net heating 
is called the gain radius, Rg, typically at ~100 km. Neutrino 
heating is provided by both accretion luminosity as the matter 
falling through the bounce shock is incorporated into the PNS 
and core neutrino luminosity as the PNS cools and deleptonizes. 

During collapse, the inner parts of the iron core collapse to 
smaller radii than the outer parts. Thus, if one begins with a 
sohd-body rotation profile at the moment the core becomes un- 
stable, as the collapse proceeds, strong differential rotation - a 
negative shear (c/fi/c/ In r, where ft is the angular velocity and 
r is the radial coordinate) profile - naturally develops (LeBlanc 
& Wilson 1970; Fryer & Heger 2000; Akiyama et al. 2003; Ott 
et al. 2003). The strong shear profile created during collapse 
persists as the shock is formed, as it stalls, and as the accretion 
phase begins. The shear continues to grow, particularly at the 
PNS surface, as material is advected inward and the whole pro- 
genitor core is accreted. This shear profile is a reservoir of free 
energy, a portion of the gravitational binding energy of collapse 
trapped as differential rotation. 

Heger et al. (2000) have computed the evolution of rotating 
stars with zero-age-main sequence masses greater than 10 
and with initial equatorial surface angular velocities up to ap- 
proximately 70% of breakup. Although spherical, these models 
account for angular momentum transport via a variety of insta- 
bilities and include the effects of the centrifugal force up to cen- 
tral neon burning. Heger et al. (2000) neglect rotational effects 
in the momentum equation and energy equation during the evo- 
lution from neon burning to core collapse. The specific angular 
momentum they predict in the roughly solid-body iron core that 
forms prior to core collapse is j ~ 10'^- 10'^ cm^ s~', with cen- 
tral angular velocities ~ 1 - 10 rad s~^ On collapse, the PNS 
is formed with a mass ^1.2 Mq and a radius of ^ 8 x 10^ cm. 
The ratio of pre-collapse to post-collapse radius for a mass el- 
ement which comprises the PNS is typically 10-20, implying, 
via angular momentum conservation, an increase in the ratio 
of pre-collapse to post-collapse angular velocity of 100-400. 
Taking intermediate values for the initial angular velocity of the 
iron core, the total energy in shear can then be estimated: 

corresponding to a spin period near the PNS surface of ~ 6 ms. 
The total energy in shear in the gain region, exterior to the PNS, 
may be less by an order of magnitude, but is replenished on a 
freefall timescale as more of the progenitor falls through the 
standoff supernova shock. Although shear energy in the range 
of 10^' - 10^^ erg is small on the scale of the binding energy 
of the resulting fully contracted neutron star and the energy 
liberated in neutrinos, it is large on the scale of the canonical 
asymptotic supernova energy and (not coincidentally) the bind- 
ing energy of the matter exterior to the gain radius, ^ 10^' erg. 
If this free energy source can be tapped on a timescale compa- 
rable to, for example, the neutrino heating timescale in the gain 
region, we expect large modifications to the dynamics. 

Any local viscous process acting in a region of differential 
rotation will transport angular momentum and dissipate the en- 
ergy stored in shear on a viscous timescale. It is this mecha- 
nism for transferring gravitational binding energy, stored dur- 
ing collapse as shear energy, to the matter in the gain region 
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that we explore in this paper. To this end, we consider sev- 
eral potential mechanisms for viscosity: microscopic shear neu- 
trino viscosity and the turbulent viscosity caused by convective 
or magnetic stresses. Generation of these turbulent stresses is 
caused by unstable entropy, composition, or angular velocity 
gradients. As an expedient, we distinguish between magnetic 
stresses caused by the magnetorotational instability (MRI) and 
the turbulent viscosity of purely hydrodynamical convection 
(see §4.2 and §4.3). We compare these agents of viscosity and 
find that the MRI and hydrodynamical convection dominate mi- 
croscopic neutrino viscosity and that the latter is hkely not ca- 
pable of modifying the dynamics of supemovae. We find that 
for fast enough initial rotation periods, dissipation of shear en- 
ergy can affect the subsequent dynamics significantly; a stalled 
shock can be revived and driven to infinity; a dud is transformed 
into an explosion. We emphasize that for the success of this 
mechanism we do not require large seed magnetic fields in the 
iron core prior to collapse. The large magnetic fields generated 
from even a vanishingly small initial field are a necessary con- 
sequence of the MRI in rapidly rotating collapse (Akiyama et 
al. 2003). 

The calculations presented in this paper are one-dimensional. 
Spherical symmetry is assumed here, but demonstrably bro- 
ken in a rotating system. We have attempted to partially mit- 
igate this issue by restricting our calculations to those in which, 
throughout the evolution, Keplerian rotational velocities are not 
attained. That is, for every model presented, at every time and 
at all radii, r2(r) < = \jGM/r'^. By imposing this restric- 
tion, the momentum equation is dominated by the gradients in 
the thermal pressure and the gravitational force; we calculate 
stars, not disks. That said, for Zishear to be in an interesting 
energy range (say > 10^*^ erg) rotation rates for the iron core 
at collapse must be in the range of a few rad s"^ . Assuming 
strict conservation of angular momentum and that the super- 
nova leaves behind a ^ 1 .4 M© neutron star with a radius of 
~ 10 km, we find that this energy requirement implies an ini- 
tial spin period for the young neutron star of at most ~10 mil- 
liseconds. Observations of radio pulsars imply that the initial 
rotation rate of neutron stars should be in the range of many 
tens to hundreds of milliseconds (e.g. Kaspi & Helfand 2002), 
factors of hundreds to tens of thousands lower in rotational en- 
ergy. Thus, although we restrict ourselves to models that are 
slowly rotating on the scale of Ok. they are rapidly rotating 
on the scale of observed neutron stars. However, caution is 
warranted in taking this argument too seriously. Gravitational 
radiation (e.g. Ostriker & Gunn 1969), magnetic dipole radi- 
ation (e.g. Pacini 1967, 1968; Gunn & Ostriker 1969; Lind- 
blom et al. 2001; Arras et al. 2003), magneto-centrifugal winds 
(Thompson 2003; Thompson, Chang, & Quataert 2004), and 
late-time fallback (see Woosley & Heger 2003 and references 
therein) can spindown young neutron stars, complicating the in- 
ference of "initial" neutron star spin periods (Kaspi & Helfand 
2002). Nevertheless, it is possible that rapid rotation is relevant 
in only a subset of all core-collapse supernovae (e.g. those pro- 
ducing magnetars; Duncan & Thompson 1992; Thompson & 
Duncan 1993). The size of that subset is unknown. 

In addition to the various processes that may slow a neutron 
star at or just after birth, it is also possible that stellar evolution 
produces slowly rotating iron cores. Recent models of Heger et 
al. (2003) and Heger et al. (2004), based on the work of Spruit 
(2002, 2003), attempt to evolve core-collapse progenitors in- 



cluding angular momentum transport via magnetic processes. 
These effects yield significantly smaller iron core specific an- 
gular momenta (j ~ lO'"*- 10'^ cm^ s~\ O < 0.1 rad s~^), im- 
plying much lower shear energy. Even smaller rotation periods 
were found in Spruit & Phinney (1998). Hence, uncertainty in 
models of rotational stellar evolution remains high. 

Similar uncertainty surrounds models of rotating core col- 
lapse and the role of rotation and magnetic fields in the mech- 
anism of supemovae. Although pioneering, the MHD calcu- 
lation of rotational core-collapse by Leblanc & Wilson (1970) 
is now largely obsolete because the neutrino microphysics, nu- 
clear equation of state, and progenitor configuration they em- 
ployed have been superseded. The jet-hke explosion their cal- 
culation produced was reconsidered with more modem input 
physics by Symbalisty (1984), who showed that only with ini- 
tial iron core spin periods much less than one second and a 
very large fossil seed field in the progenitor iron core, could 
jet-driven explosions be obtained. 

The recent work of Akiyama et al. (2003) considered the po- 
tential effects of the MRI in stellar collapse. They showed that 
much of the post-bounce core collapse environment is unstable 
to the MRI and that large saturation magnetic fields are natu- 
rally produced. Energy deposition via local turbulent dissipa- 
tion was neglected in their study. Instead, Akiyama et al. (2003) 
argued that the very large field strengths generated by the MRI 
create a jet with sufficient MHD luminosity to explode the star 
as in the models of Khokhlov et al. (1999). This argument as- 
sumes that the magnetic field generated by the MRI can form 
the organized large-scale fields required for coUimation and jet 
formation. A limitation of the models presented by Akiyama 
et al. (2003) is that they were not self-consistent. They did not 
include rotation in the dynamics and because they neglected an- 
gular momentum transport, the energy available for the MHD 
jets they consider may have been overestimated. 

Fryer & Heger (2000) considered rotating core-collapse in 
two spatial dimensions, and recently Fryer & Warren (2004) 
extended these rotating models to three dimensions. The only 
effect of rotation on the explosion mechanism that Fryer & 
Heger (2000) and Fryer & Warren (2004) discuss is the sup- 
pression of convection in the equatorial plane due to the sta- 
bilizing centrifugal force. However, for stabilization to occur, 
rotation must be important enough that if the calculation had 
been performed in MHD the equatorial region would have been 
destabilized by the MRI. Thus, the stabiUzation that occurred 
in their hydrodynamic calculations is somewhat artificial and 
it is a signal that the simulations must be performed in MHD. 
All models computed in Fryer & Heger (2000) and Fryer & 
Warren (2004) explode. For this reason, Fryer & Heger (2000) 
and Fryer & Warren (2004) did not assess the role of viscous 
heating, the competition between different possible sources of 
viscosity, or the potential role of rotation itself in generating 
explosions. Finally, Fryer & Warren (2004) concluded that the 
saturation magnetic field strength throughout the post-collapse 
structure is orders of magnitude lower than the field estimated 
in Akiyama et al. (2003). We address this conclusion in §5 and 
find that large fields in the range 10^"*- 10'* G can be typical in 
rotating core collapse ( Duncan & Thompson 1992; Thompson 
& Duncan 1993; in agreement with Akiyama et al. 2003). 

3. THE MODEL 
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The core-collapse supernova code we have developed and 
that we employ here, SESAME^, is described in detail in 
Burrows et al. (2000) and Thompson et al. (2003). The hy- 
drodynamics scheme is Lagrangean, Newtonian, explicit, and 
employs artificial viscosity for shock resolution. The radia- 
tion transport algorithm solves the order V/c Boltzmann equa- 
tion (Eastman & Pinto 1993) for three neutrino species: Ve, 
Ue, and 'V^", employing the standard approximation that the 
non-electron neutrino species can be treated collectively. It 
is fully spectral, employs the tangent-ray technique for angu- 
lar discretization of the radiation field, and is based on the 
Feautrier variables. The neutrino physics employed is stan- 
dard and described in Bruenn (1985), Thompson (2002), Bur- 
rows & Thompson (2002), and Thompson et al. (2003). No- 
table additions include nucleon-nucleon bremsstrahlung (Bur- 
rows et al. 2000; Thompson et al. 2000), and weak-magnetism 
and recoil corrections to both the charged-current neutrino- 
nucleon absorption/emission processes and elastic neutrino- 
nucleon scattering (Horowitz 1997; Horowitz 2002). Inelas- 
tic neutrino-electron/positron scattering is calculated with a 
fast and accurate explicit scheme described in Thompson et 
al. (2003). The equations of radiation hydrodynamics are cou- 
pled to an efficient tabular version of the Lattimer-Swesty high- 
density nuclear equation of state (for p > lO^gcm"^) and the 
Helmholtz equation of state (for p < lO^gcm"^; Timmes & 
Amett 1999; Timmes & Swesty 2000; Thompson et al. 2003).* 
AH of the models presented here are calculated with 20 energy 
groups per neutrino species and 500 Lagrangean mass zones. 

3.1. Additions to the Standard Algorithm 

To the standard set of Newtonian equations for Lagrangean 
hydrodynamics in one spatial dimension, we add rotation via 
the term fej^ /r' in the momentum equation, where /' is the 
specific angular momentum, r is the radial coordinate, and 
/©(= 2/3) results from taking the angular average for one di- 
mensional parameterizations of rotation (Heger et al. 2000). 

Angular momentum transport is calculated by solving the 
equation 



The dissipation of shear energy is handled by adding a term 
to our Lagrangean equation for the specific internal energy (e); 
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(3) 



where the actual form of ^ and the resulting viscous energy 
deposition rate (q^) depends on the specific viscous process 
considered. Note that we dissipate shear energy locally. For 
a source of microscopic viscosity this is appropriate. However, 
for a turbulent viscosity like the MRI, it may be that energy 
stored in the turbulent magnetic field can be transported and de- 
posited non-locally. This would relieve the requirement in the 
models presented here that a region of high local viscous dis- 
sipation must correspond to a region of large shear. However, 
a multi-dimensional MHD treatment is required to address this 
question fuUy in the supernova context. 

3.2. The Progenitors 

We consider two progenitors. The first is the non-rotating 
1 1 Mq progenitor of Woosley & Weaver (1995). This model is 
zoned out to 1.45 Mq, corresponding to a radius of ~ 4000 km. 
We impose rotation on the core of the progenitor, taking 

where r is the spherical radial coordinate and not the cylindrical 
radial coordinate. Equation (4) implies that the rotation profile 
is roughly sohd-body out to Rq. Comparing with the rotat- 
ing stellar progenitors of Heger et al. (2000), eq. (4) provides 
a reasonable approximation to rt(r) with Rq near ~ 1000 km. 
Our models span spin periods in the range 1 ^ Pq ^ 10 sec- 
onds {Pq - In/rio). As discussed in §2, we restrict ourselves to 
models in which Q{r) < Qk at all times and at all radii during 
a given calculation. For our 1 1 Mq model and for f2(r) set by 
eq. (4) with Rq = 1000 km, we find that Pq > 1 .25 s satisfies this 
constraint. This sets our lower bound on Pq. The upper bound 
on Pq is set by the fact that our spherical models with Pq > 8 
s and Rq = 1000 km are nearly indistinguishable from models 
with Po = oo. In all models with this 1 1 Mq progenitor we take 
Rq = 1000 km. 

The second progenitor we consider is the rotating 15 M© 
model E15A of Heger et al. (2000), which we zone out to an 
enclosed mass of 1.65 Mq. With this model we can compare 
directly to the 2D axisyrmnetric core-collapse calculations of 
Fryer & Heger (2000) and the 3D calculations of Fryer & War- 
ren (2004). 

Figure 1 shows J7(r) for the progenitor models considered 
here. Solid lines are for the 11 Mq progenitor with Pq = 1.25, 
2, 3, 4, 5, and 8 seconds. The dashed line shows il(r) for the 
15 Mq model from Heger et al. (2000). Throughout this paper 
we use "E15A" to refer to the 15 M© progenitor of Heger et 
al. (2000) and "Pq =" in referring to the 11 M© models. 

^ Spherical Explicit/Implicit Supernova Algorithm for Multi-Group/Multi-Angle Explosion Simulations. 

* A recent comparison between the Lattimer-Swesty (LS)EOS in NSE at moderate densities and temperatures (characteristic of the gain region in mod- 
els of supemovae) with other NSE EOSs indicates a discrepancy of potential importance. Free nucleon and a-particle nvmiber fractions as well as en- 
tropy profiles disagree at the ~ 10 - 30% level in some regions of thermodynamic space. The comparison is documented by F. X. Timmes online at 
http://www.cococubed.com/code _pages/eos_supemova.shtml. Janka et al. (2004) have compared supernova models computed with the LSEOS with otherwise iden- 
tical models computed with the high-density EOS of Shen et al. (1998). They find quantitative differences between the two models, but they do not find that those 
differences lead to quahtatively different model outcomes; both models fail to yield explosions. Because these models do not show large qualitative differences and 
because in this paper we make a relative comparison between models without rotation, with rotation, and with viscosity and angular momentum transport, we have 
employed the standard LSEOS. 

^ Our use of to denote "viscosity" alleviates confusion with 'V", which we use to denote neutrinos or neutrino quantities. 
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pr^ dr 



drj' 
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where D/Dt = 9/9f + V-V, p is the mass density, is the 
angular velocity, and ^ is the shear viscosity.^ Equation (2) 
is solved explicitly in operator-split fashion from the rest of 
the radiation-hydrodynamical algorithm. The calculations pre- 
sented in this paper conserve total angular momentum at the 
~ 0.05% level; maximum fractional deviations in the total an- 
gular momentum do not exceed ~ 0.15%. During a given cal- 
culation, depending on the magnitude of the viscosity, the diffu- 
sive timestep (At < min[(Ar)^/2^]) becomes shorter than the 
Courant-Umited hydrodynamical timestep and we sub-cycle the 
calculation of the update to j. 
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4. VISCOSITY 

In this section, we discuss the pure microscopic shear vis- 
cosity caused by neutrinos and two forms of turbulent viscos- 
ity: the magnetorotational instability and hydrodynamic con- 
vection. We compare the magnitude of these viscous processes 
in §6 in the context of the rotating models of core collapse pre- 
sented in §5. In §7, our results for rotating models including 
viscous dissipation are discussed. 

4.1. Neutrino Viscosity 

When neutrinos are fully diffusive, the shear viscosity is 
given by (van den Horn & van Weert 1984; Burrows & Lat- 
timer 1988; Thompson & Duncan 1993) 



where n is the epicyclic frequency. 



(10) 



<ilnr' 

and is the Brunt- Vaisala frequency. In the supernova context, 
material may be driven convectively unstable by negative gra- 
dients in entropy (s) or in lepton number (F/). For our purposes 
here, it is sufficient to replace Y[ with the electron fraction (Yg), 
thereby neglecting the electron neutrino fraction in regions 
where this species is trapped and diffusive. The expression for 
the Brunt- Vaisala frequency used in this paper is then (Lattimer 
& Mazurek 1981) 



ds 



P,Ye 



ds 1 dP 
dr^P We 




^ , (11) 
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where E^, is the neutrino energy density, {X^) is the average 
neutrino mean-free path. 



(A,.) = { deu£lX^{eu)Ju{£u) \ I { i de„elJ^(e^) 
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diifv{li,e„), 



(6) 



(7) 



where g is the effective gravity and 7 = d\nP /d\np\s. When 
magnetic fields are included in the stability analysis of differ- 
entially rotating bodies, a "generalized" Solberg-H0iland crite- 
rion is obtained (Balbus & Hawley 1991; 1992a,b; 1994; 1998). 
The condition 

A^'+Ti— <0, (12) 
fl Inr 

replaces eq. (9) as the criterion for instability in the equato- 
rial region. Ignoring gradients in entropy and composition, a 
magnetized differentially rotating fluid is unstable to the mag- 
netorotational instability if dff'/dlnr < 0. The linear growth 
rate for the fastest growing mode is (Balbus & Hawley 1991) 

1 dn 

i MRI = 



dlnr 



(13) 



To rough, but useful, approximation, Fmri ~ f2. In a sub-sonic, 
steady-state accretion flow with net heating, as in the gain re- 
gion just interior to the shock in the supernova context, the 
entropy gradient is negative, and compositional gradients are 
insufficient to stabilize this region. Because the process of col- 
lapse, along with angular momentum conservation implies neg- 
ative shear in the gain region, eq. (12) implies that as long as 
Ti/ < Fmri this region must be unstable to the MRI. 

Our form for the viscosity due to the MRI is motivated by 
consideration of the MaxweU stress (trd,); 



tr4> = BrB^/A-K ~ B^/Att = VIp, 



(14) 



fv is the invariant neutrino phase-space distribution function, 
£u is the neutrino energy, and fi = cos9. 

Such a prescription for the viscous effects of neutrinos ap- 
plies only in the region of diffusive neutrino transport, below 
the neutrinospheres, where (A^) <C r. In the optically thin re- 
gion, as (A,y) 00, eq. (5) is inapplicable. Nearest the neutri- 
nospheres, will be large, but outside this region the viscosity 
is not well defined. However, the effect of the streaming neutri- 
nos interacting with the background differential rotation profile 
and the distribution of hydrodynamic (convective) and magne- 
tohydrodynamic (magneto-convective) fluctuations can still be 
quantified. Neutrinos propagating through the cooling and gain 
regions damp fluctuations with a characteristic rate (Jedamzik 
et al. 1998; Agol & Krolik 1998) 

r.--^. (8) 

For example, if a convective mode of wave number k in the gain 

region has linear growth rate given \N^ \ , then if > \/\N^\ 
the convective mode will be damped. Similarly, if the MRI 
grows in the gain region on a timescale ~ if > the 
linear MRI mode will also be damped. 

4.2. The Magnetorotational Instability 

In the equatorial region, neglecting gradients in the 0- 
direction^, the Solberg-H0iland condition for instability in the 
presence of rotation, but without magnetic fields is given by 
(Tassoul 1978) 

N^ + K^<0, (9) 

* The Solberg-H0iland criterion for instability contains latitudinal gradients that for a star may be as large as those in the radial direction. For the purposes of the 
discussion here, and in the spirit of the one-dimensional calculations presented in this paper, we retain only the latter. 

' Akiyama et al. (2003) estimate the MRI viscous timescale as ryisc ~ B^Pu/C^Br^l), where Pm = P/(B^/87r) (see their eq. 32). Because /3m S> 1 in a star, they 
neglect the redistribution of angular momentum by magnetic stresses. This is incorrect. In estimating the a parameter, they take a ~ BrB^/(4irP), and it is in this 
way that /3m enters their Tvisc- Their estimate for a, however, is relevant for disks only. In the stellar context, tVisc ~ (o^^T^ir/H)^. For most of our models this 
estimate yields tvisc ~ 0.1 s at Ri/ and tvisc ~ 1 s at Kg and so angular momentum redistribution cannot be neglected. 



where Va = B/^4np is the Alfven velocity. The viscous shear 
stress is given by ^p\dfl/d]nr\ ~ ^p^. Equating the MaxweU 
stress and the shear stress we find that ^mri ^ ^a/^- the 
MRI, VaA: ^ n for the fastest growing mode of wave number k. 
Taking A: ~ 1 /H, where H is the pressure scale height, and sub- 
stituting, we find that ^mri ~ i/^il, up to factors ofdla^l/dlnr. 
We introduce a constant of order unity, an a parameter (typi- 
cally 0.1), and take 

^Mm = aH^n. (15) 

Equation (15) is appropriate for the form of the viscous stress 
regardless of the ratio V^/c^.^ In all models considered in this 
work, is less than or much less than Cj and Vkcp- This con- 
dition emphasizes the fact that we consider stars, not disks. In 
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accretion disks, H = Cj/il. Substituting into eq. (15) one ob- 
tains the classic form for the shear stress in accretion disks from 
Shakura & Sunyaev (1973), ^ = acsH. Another way to see that 
eq. (15) is the appropriate form for the viscosity is to recognize 
that any turbulent viscosity can be estimated as the product of a 
turbulent velocity (Vturb) and a correlation length (L) for fluctu- 
ations. In the case of a disk the adiabatic sound speed c, is Vturb 
and L is the pressure scale height, H. In a star, with HVl less 
than or much less than Cj, the shear across a scale height HQ. 
takes the place of Cs as Vturb- The correlation length for the MRI 
remains the pressure scale height and eq. (15) follows. 

Beginning with an arbitrarily small seed field, the MRI in- 
creases the magnetic field strength exponentially on a timescale 
^ f2"^ In models of accretion disks, the MRI in ideal MHD 
saturates such that the toroidal component of the Alfven speed 
comes into rough equipartition with the sound speed. In the 
stellar context considered here, the Alfven speed should satu- 
rate at of order the rotation speed, 

Bsat, ~ (47rp)i/2y^. (16) 
The value of the toroidal saturation magnetic energy density 
(^Sat^/^"") ^® somewhat less (a factor of ~10) than the 
value impUed by eq. (16) (see, e.g., Hawley et al. 1996 and 
Stone et al. 1996 for results in the accretion context). A re- 
duction of this magnitude in the energy density implies a rather 
small correction to the field strength Bsat^ ■ Further possible 
corrections to eq. (16) are order unity and include factors of 
d\nQ,/dh\r (Akiyama et al. 2003). The magnitude of the sat- 
uration energy density in the poloidal direction is roughly an 
order of magnitude smaller than Bsat^ (e.g. Balbus & Hawley 
1998). 

4.3. Convective Viscosity 

The Brunt- Vaisala frequency is always negative in the gain 
region. Convection, driven by neutrino heating, results (Herant 
et al. 1994; Burrows, Hayes, & Fryxell 1995; Janka & Miiller 
1996; Mezzacappa et al. 1998; Buras et al. 2003). If convection 
has time to develop fully, and the Reynolds number is high, the 
presence of turbulent convection may result in the dissipation 
of shear energy in much the same way as turbulence caused by 
the MRI. The convective viscosity (^con) can be naively esti- 
mated as the product of the turbulent velocity (in this case the 
convective velocity Vcon) and the correlation length: 

eCon~//Vcon/3. (17) 

The convective velocity can be estimated from mixing length 
theory: 

Vcon = 2gl{5p/p), (18) 

where g = -p~^dP/dr = P/ pH, I is the mixing length, and 
6p = {dp/dr- /S.p/l)l. In the gain region, Vcon is typically 
^ 5 X 10^ — 10^ cm s^' . Unfortunately, this estimate is very 
rough because even the sign of the Reynolds stress, which dic- 
tates the direction of angular momentum transport in a con- 
vective region, is uncertain and may depend sensitively on the 
Rossby number, the ratio of the rotation period to the convective 
turnover time (Chan 2001; Kapyla, Korpi, & Tuominen 2003). 
In general, for Rossby numbers less than unity the Reynolds 
stress is positive and angular momentum transport is outward. 
However, this statement depends on latitude and the magnitude 
of the heat flux through the convection zone. In the supernova 
context, particularly for rotation rates with V^ <C Cj in the gain 



region (and deep in the PNS interior), we expect the Rossby 
number to be greater than unity (Thompson & Duncan 1993). 

It is not strictly correct to speak of two different mechanisms 
for turbulent convective viscosity; the MRI and hydrodynami- 
cal convection come together. There is just one turbulent stress 
tensor, containing Reynolds and Maxwell correlations, of the 
form {urU^ - ua,ua^) , where u is the fluctuation velocity and 
Ma, = Bj/ sJ\-Kp is the fluctuation Alfven speed in direction i\ 
and convection itself will, of course, amplify magnetic fields 
even in the absence of the MRI, thereby generating apprecia- 
ble magnetic stresses (Duncan & Thompson 1992; Thompson 
& Duncan 1993). But one term in the stress tensor may domi- 
nate the other and for the purposes of this paper we distinguish 
between purely hydrodynamical convection (Fcon 3> Fmri) and 
the case where the Maxwell stress is comparable to or domi- 
nates the Reynolds stress (see Narayan et al. 2002 for a more 
general discussion). 

5. CORE COLLAPSE WITH ROTATION, BUT WITHOUT 
VISCOUS DISSIPATION 

In this section we describe results from a set of simulations 
wifli the 11 Mq progenitor wifli = 1-25, 2, 3, 4, 5, and 8 
seconds as well as results with model E15A from Heger et 
al. (2000) (see §3.2). In these models we include the centrifugal 
force, but neglect diffusive angular momentum transport. We 
do this in order to separate the effects of the centrifugal force 
from the effects of viscous dissipation. Table 1 summarizes 
a few properties of these rotating models. We list the model 
names, the shear energy at bounce and 500 ms after bounce, 
and the ratio of total rotational energy to total gravitational en- 
ergy (/3^°') initially, at bounce, and 500 ms after bounce. In §7, 
we describe results from these models when viscous dissipation 
is included. 

Figure 2 shows the neutrino luminosity (left panels, in units 
of 10^^ ergs"') and average rms neutrino energy (right panels, 
in MeV) at infinity as a function of time relative to hydrody- 
namical bounce for each of our rotating models. SoUd Unes are 
for the 1 1 M0 progenitor of Woosley & Weaver (1995). Dashed 
lines are for model E15A from Heger et al. (2000). Specifica- 
tions for the models are given in §3 and Table 1. Note that in 
none of the models shown here is a centrifugal barrier reached 
during collapse. All models undergo bounce near the nuclear 
phase transition in the inner core and have sub-Keplerian ro- 
tational velocities throughout the computational domain. For 
clarity of presentation, the peak of the breakout pulse of elec- 
tron neutrinos is not shown in the upper-left panel. The most 
rapid rotator (fo = 1-25 s) yields Ly^'* ~ 3.1 x 10^^ erg s"', 
whereas both E15A and the model with Pq = 2 s have I^f^ ~ 
2.6 X 10^3 erg s'^ . AU other models have I^f^ ~ 2.4 x 10^^ erg 
s"'. 

As expected, we find that the higher the initial rotation rate, 
the lower the core temperatures of the PNS formed just after 
bounce. On average, after electron neutrino breakout, this effect 
produces lower core and {ej) for shorter Pq. This reflects 
the sensitivity of the neutrino emissivities to the local temper- 
ature. The fractional differences in (ei,) between the model 
with Pq = 1 -25 s and a non-rotating model 200 ms after bounce 
are approximately 15%, 17%, and 30% for v^, Dg, and Vg, re- 
spectively. The same comparison for L^, yields fractional dif- 
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ferences of 75%, 63%, and 33% for L^^, Lp^, and L^,^, respec- 
tively. The difference in L^^ at 100 ms after bounce between 
our slowest and fastest rotators is a factor of ^ 6. The neutrino 
spectral characteristics change significantly and systematically 
as a function of Pq. Despite the fact that model E15A has rapid 
rotation, it has a much larger L^^ after breakout than, for exam- 
ple, the Po = 2 s 11 Mq model, even though they have similar 
initial f2(r) profiles (see Fig. 1). This is due to the extended 
density profile of model E15A relative to the 1 1 Mq model and 
the associated larger accretion luminosity. 

Figure 3 shows radial profiles of a number of quantities in 
a subset of the models of Fig. 2 at a time-slice 105 ms after 
bounce. The upper left panel shows the mass accretion rate for 
the models with Po = 1.25 (dotted), 2 (dashed), 3 (dot-dashed), 
and 8 seconds (solid). The shock is clearly visible in each pro- 
file, the more rapid rotators having larger shock radii 105 ms 
after bounce. Higher negative M imphes larger PNS accretion 
rate and a larger accretion luminosity. This, coupled with the 
fact that the faster rotators achieve lower core temperatures at 
bounce, explains why the total luminosity is smaller for smaller 
Pq in Fig. 2. The upper right panel shows fl (positive numbers, 
in units of 10^ rad s"^) and d^/dlnr (mostly negative numbers, 
in units of 10^ rad s"'). The same line styles are employed. 
Note that the surface of the PNS (the neutrinosphere at the peak 
of the i^e spectrum) in this epoch lies at approximately 50 - 70 
km (depending upon Po). Thus, if the mantle were ejected at 
this instant, assuming strict angular momentum conservation, 
all models with Pq ^ 8 s form 10 km neutron stars with P < 10 
ms. In fact, the model with Pq = 8 s would have a final spin 
period of ~ 7 ms. The model with Pq = 2 s would have a final 
angular velocity over lO'' rad s~^ 

The lower left panel of Fig. 3 shows the saturation magnetic 
field strength if the MRI operates everywhere, obtained simply 
by equating the magnetic energy density with the azimuthal ki- 
netic energy density (as in eq. 16). Figure 3 shows that at radii 
near 100 km magnetic fields in the range ~ lO''*- 10^^ G 
would be generic to relatively slowly rotating progenitors (here 
typified by the Pq = 8 s calculation, solid line). For rapid ro- 
tators, much of the core has field sfi-engths of ~ 10^* G and 
above. The results presented here for Bsat^ are in rough agree- 
ment with the models of Akiyama et al. (2003), but much higher 
than those presented in Fryer & Warren (2004). Comparing our 
Fig. 3 with their Fig. 13 we find that our field strengths are gen- 
erally larger by 4-5 orders of magnitude. 

Although we have not included viscous dissipation in any of 
the models presented in this section, in the lower right panel of 
Fig. 3 we plot the magnitude of the viscous dissipation rate, as- 
suming that the MRI operates everywhere. We combine eqs. (3) 
and (15) to obtain 

.MHi = ^M.(^) =a//^0(^j , (19) 

where we take a = 0.1. The solid box shows a range of rep- 
resentative net neutrino energy deposition rates in the region 
exterior to the PNS (60 < r < 200 km). Typical net neutrino 
heating rates in the gain region are ~ 100 MeV nucleon"^ 
s"'~9.6 X 10'^ erg g"' s"'. The box brackets the range lO'^ ^ < 
^MRi < 10^°^ erg g"' s"' . Of course, between the surface of the 
PNS and the shock the net heating rate varies considerably, go- 
ing from negative to positive at the gain radius Rg. This panel is 
meant merely to underscore that a first comparison shows that 



the neutrino heating rate and the viscous heating rate are com- 
parable for rapid rotators. If enough shear is present so that the 
energy deposition rate imphed by the lower right panel of Fig. 3 
can be maintained, the dynamics of the collapsed core may be 
significantly affected (§5.1). This is particularly true for more 
rapid rotation because the net neutrino heating rate in the gain 
region decreases as Cl increases and the viscous heating rate in- 
creases as fP. Thus, the relative importance of viscous heating 
with respect to neutrino heating is compounded as the rotation 
rate increases. 

For clarity of presentation we have not shown results for 
model E15A in Fig. 3. The saturation magnetic field strength 
and viscous dissipation rate for this model are very similar to 
the model with Pq = 2 s in the lower left and right panels, re- 
spectively. Model E15A has a larger negative M at P^ and a 
larger R,j than model Pq = 2 s, but similar shock radius. Its shear 
and angular velocity profiles closely ttack those of model Pq = 2 
s for r > 20 km. For r < 20km model E15A has O > lO'' rad 
s"' , greater than the corresponding number for model Po = 1 .25 
s. These differences between model E15A and model Po = 2 s 
reflect the difference in initial il{r) for r < 600 km (Fig. 1) and 
the fact that the density profile for model E15A is much broader 
than that for Pq = 2 s. 

5.1. Why Supernovae do not Explode in ID 

Figure 3 shows that the viscous heating via the MRI can be 
comparable to the net neuttino heating in the gain region. In 
§7, we will show that this effect is sufficient to yield explo- 
sions in some models. In order to address the question of why 
these models succeed, it is instructive first to consider why non- 
rotating spherical models with the very best physics fail. An 
elucidating and physical discussion of the stalled shock dynam- 
ics can be found in Janka (2001). We do not attempt to repro- 
duce the detailed arguments presented there, but instead offer 
a semi-quantitative quasi-global criterion for explosion in the 
spirit of the arguments presented in Thompson (2000). 

First, we define a radial advective timescale across a pres- 
sure scale height TAdv = H/Vr and a timescale for net neuttino 
heating r^^ = {P/ p)/qv, where P is the thermal pressure, p is 
the mass density, q,^ = //^ -C^, and H,^ and Ci, are the spe- 
cific neutrino heating and cooling rates, respectively. Roughly 
speaking, steady-state accretion is maintained (explosions do 
not develop) in simulations of core-collapse in one dimension 
because TAdv < in the gain region. If this were not the case 
(that is, if TAdv > Tq^^), the thermal pressure would increase be- 
tween Pg and Psh (or, not decrease as quickly as material is 
advected into the coohng region) and the shock would move 
outward, finding a new equilibrium. This condition, TAdv > Tqu > 
is a necessary but not sufficient condition to guarantee explo- 
sion. One must maintain the condition TAdv > for a time 
sufficient to deposit roughly the binding energy of the overly- 
ing mantle, 10^' erg. Thus, if the condition TAdv > is met 
for ti<t < tf and 

E^^s< J" J<}^^pd'rdt, (20) 

then the mantle is unbound and an explosion results. The 
statement of eq. (20) is equivalent to saying that the matter 
achieves the "escape" temperature, obtained by equating the 
local specific thermal energy with the gravitational potential 
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energy (Burrows & Goshy 1993; Burrows, Hayes, & Fryxell 
1995). 

We illustrate these basic numbers in Fig. 4, which shows Tg^ 
(solid line), TAdv (dashed line), rc^ = (P/p)/Cu (dotted Une), 
and th^ = iP/p)/H^ (dot-dashed line) for a slowly rotating ini- 
tial progenitor (Pq = 8 s, see Fig. 2) 130 ms after bounce. The Vg 
neutrinosphere (R^J and the gain radius (Rg, where tc„ crosses 
th^ and r^^ diverges) are indicated by arrows at this snapshot in 
time. Note that throughout the gain region (r > Rg) TAdv < Tq^ ■ 
The condition for explosion we suggest does not obtain. Fol- 
lowing the evolution in this model for another 400 ms, Tq^ / TAdv 
is always greater than 1 at its minimum (which usually occurs 
very near, but exterior to, Rg). We find very similar absolute val- 
ues for the various timescales and ratios t^^ / TAdv for all models 
that fail to explode. Of course, one may argue that the actual 
value of, for example, TAdv, depends on its definition and that 
it may differ by order unity from that presented in Fig. 4. For 
example, we define TAdv as H/Vr rather than r/Vr because 
depends fairly sensitively on radius near Rg since the cooling 
rates are highly dependent on the local temperature. Thus, the 
distance over which q^, doubles near Rg is closer to H than r. 
This motivates our definition of TAdv. but there is no reason a 
priori that TAdv should not differ from H/Vr by, say, a factor of 
~ 2. In addition, Tg^ as we have defined it corresponds to an 
e-folding time for the local pressure, a relatively large change 
in the gain region. It may be that well before a full e-folding of 
P is accomplished the system adjusts. This is particularly true 
because the dynamical timescale (Toyn =H/cs or r/d) is much 
shorter than all other timescales. Thus, we cannot rule out or- 
der unity corrections to Tg^ . For these reasons, that the ratio 
Tq^ I TAdv is ^ 1 .3 at r ^ 125 km does not necessarily imply that 
if the net heating rate is increased by ~ 30% the shock must be- 
gin to move outward. However, the ratio r^,^ / TAdv is certainly 
indicative of the magnitude of the deficit in heating that must be 
paid for explosion. Remarkably, we show in §7 that the crite- 
rion advocated here is quantitatively useful for assessing when 
explosions will develop. 

Although simplistic, the picture presented here, which cor- 
responds in rough terms with the conditions derived in Janka 
(2001), yields predictive power In all of our spherical, non- 
rotating models (see Thompson et al. 2003) TAdv is always less 
than Tq^ in the gain region over the full post-bounce evolu- 
tion. Similarly, in model E15A and the models with fo = 2, 
3, 4, 5, and 8 seconds, despite the wide range of mass accre- 
tion rates and post-bounce structures, the shock stalls in each 
model for the duration of the calculation and T^iw/Tq^ is al- 
ways less than unity (when dissipation is not included).^" Ex- 
ceptions to this rule occur only during very dynamic episodes. 
For example, the instant the Si shell encounters the shock, the 
density decreases, M decreases, and the shock moves out. In 
the instant between the two epochs of quasi- steady- state accre- 
tion, TAdv can be greater than t^^, but this is a transient and 
although the shock moves outward in radius, no explosion is 
obtained. In marked contrast, if one injects energy into the gain 
region, say via viscous dissipation with rate ^'mri, then the rele- 
vant comparison is between TAdv and Tg^^^ = {P / p)/ qioi, where 
^TOT = Hi,-Cu+qMm- Adding ^mri to qi, decreases the net heat- 
ing timescale from t^^ to t^tot- ^ ^mri is large enough to yield 

The model with Pq = 1-25 s never develops a well-defined gain region and the shock never stalls completely. The centrifugal support (the decrease in the grav- 
itational potential) as the shock is launched is sufficient to keep the shock moving, even though Keplerian rotational velocities are not obtained anywhere in the 
calculation. Matter velocities never exceed 2 x 10^ cm s~' during the calculation and the internal energy of the ejected matter is small because there is little, if any, 
neutrino heating as the explosion develops. Roughly estimating the asymptotic energy, we obtain < 10^ erg. 



'^?TOr < "^Adv and condition eq. (20) is met, an explosion results. 
We show in §7 that this mechanism succeeds. 



6. COMPARING VISCOUS PROCESSES 

In this section we compare quantitatively the magnitude of 

the neutrino viscosity with the turbulent viscosity of the MRI 
in the dense protoneutron star core where neutrinos are diffu- 
sive. We also compare the linear growth rate for purely hy- 
drodynamical convection and the MRI in the semi-transparent 
gain region with the neutrino damping rate. There have been 
a number of papers concerning the possibility that convection 
persists in PNS interiors during the accretion/explosion phase 
and afterwards, during Kelvin-Helmholtz cooling (Mayle 1985; 
Mayle & Wilson 1988; Wilson & Mayle 1988; Burrows 1987; 
Burrows & Lattimer 1988; Thompson & Duncan 1993; Keil & 
Janka 1995; Pons et al. 1999; Thompson & Murray 2001). If 
the neutrino viscosity (^^) is comparable to or larger than the 
viscosity of the MRI (^mri) or convection, any turbulent fluc- 
tuations generated by these instabilities may be damped. Sim- 
ilarly, convective or magneto-convective modes in the gain re- 
gion may also be damped if the neutrino damping rate (r,^) is 
larger than the linear growth rates for convection {\/\N^\) or the 
MRI (Tmri ~ O). 

In Fig. 5 we plot ^mri, ^v, Tmri, and \/\N^ \ for the model 
with Po = 2 s at a time 105 ms after bounce, as in Fig. 3. The 
left panel shows for each neutrino species and ^mri in the 
inner region of the PNS, inside the neutrinospheres for and 
Ve neutrinos at their respective spectral peaks. In this diffusive 
regime (Aj,) is less than or much less than H and a compari- 
son between the neutrino viscosity and the turbulent viscosity 
of the MRI is appropriate. For the purposes of comparison we 
have taken ^mri = a/Z^Q with a = 0.1. The viscosity caused by 
the MRI clearly dominates neutrino viscosity by 2-3 orders of 
magnitude over the entire profile. Thus, within the PNS, when 
the condition for instability to the MRI is met (eq. 12), the MRI 
will operate and dominate the viscosity. Even for the slowest 
rotator considered here (Pq = 8 s), ^mri is greater than by an 
order of magnitude or more throughout the PNS interior. 

In the right panel of Fig. 5 we plot Ty for each neutrino 
species and Fmri as a function of radius in the semi-transparent 
region from roughly the surface of the PNS out to the shock 
radius. Again, the linear growth rate of the MRI dominates the 
neutrino damping rates. At r ~ 50 km the ratio Tmbi/^v — 10. 
In the gain region (r ^ 150 km), where the comparison is 
most important (that is where condition eq. 12 is always met) 
rMRi/r^ ~ 10^ For the model with Pq = 8 s, Fmri/F^ ~ 100 
in the gain region. For comparison, we also show the linear 
growth rate for convection {yJ\lsp-\, dot-dashed line). As Pq 
decreases in our progenitor, \/\N^\ in the gain region stays 
roughly constant while Fmri decreases. In the model with 
Po = 8 s, x/Pv^I/Fmri ~ 10 at f ~ 100 ms after bounce, at 
r ~ 150 km. Importantly, both Fmri and i/|A^| dominate r,^ 
and we therefore expect well-developed, high Reynolds num- 
ber convection/magneto-convection throughout the gain region 
in multi-dimensional models. 
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Although the profiles evolve considerably in time, in gen- 
eral, we find that much of the PNS core and all of the gain 
region is unstable to the MRl. In the cooling region just above 
the PNS (50 km <r< 100 km), A^^ in eq. (12) is sufficiently 
large and positive to counter much of the negative shear profile 
and the MRI is stabilized. However, we expect that a multi- 
dimensional simulation of the MRI in this context may show 
magneto-convective plumes and filaments penetrating the cool- 
ing region, as in the purely hydrodynamical calculations of 
Burrows, Hayes, & Fryxell (1995). Deeper in the PNS, be- 
neath the neutrinosphere, the profile is unstable. The region 
20 km ^r<50 km is unstable to the MRI throughout most of 
the post-bounce epoch in the Pq = 2 s model. For slower ini- 
tial rotation, this region of instability shrinks in radial extent. 
Imposing n(r) as in eq. (4), we find that the very deep inte- 
rior (f < 10 km) has positive or near- zero shear (see upper right 
panel of Fig. 3) and is stable to the MRI. 

The region of very strong negative shear between 10 < r < 20 
km in the upper right panel of Fig. 3 contains strong stabilizing 
entropy gradients and is stable over most of the evolution in all 
of the P() models considered. Importantly, even though eq. ( 1 2) 
indicates that this region is stable to the MRI and convection, 
it is still possible to transport angular momentum and dissipate 
shear energy, albeit on a modified timescale (Spruit 2002). In 
stably stratified regions with differential rotation and a thermal 
diffusivity, k, Spruit (2002) argues that magnetic interchange 
instabihties lead to an effective viscosity of (see his eq. 32) 

where we have assumed a correlation length of H instead of 
r. Estimating with k = c{\^) /3 and comparing with the left 
panel of Fig. 5, we find that in the region 10 < r < 20 km 
is less than ^mri by approximately one order of magnitude and 
greater than by approximately two orders of magnitude. This 
implies a viscous timescale of r^^ = r^/^s ~ 1 - 5 seconds dur- 
ing the post-bounce epoch considered here. Interestingly, t^^ 
may be different (possibly longer) than the Kelvin-Heknholtz 
timescale for late-time cooUng. 

Although the comparison of Fig. 5 shows that the two parts of 
the generalized Solberg-H0iland criterion are of similar magni- 
tude, it is difficult to estimate the relative importance of convec- 
tion versus the MRI in the non-linear regime without recourse 
to multi-dimensional simulations. Taking our estimates of the 
turbulent viscosity for each process (eqs. 15 and 17) we have 

^Con _ 1 / Vcon \ ^22^ 

Taking Vcon ~ 5 x 10^ cms"' as representative of the results 
from multi-dimensional simulations and HQ ~ lO' cm s~' 
(P() < 3 s) we see that for a = 0.1, ^con and ^mri are of the same 
order of magnitude. We reiterate, however, that even the sign 
of the convective stress is uncertain (§4.3) so the comparison 
presented here is quite simplistic. 

7. CORE COLLAPSE WITH ROTATION & VISCOUS 
DISSIPATION 

We now provide results from a subset of the rotating mod- 
els considered in §5, including viscous dissipation and angular 
momentum transport. 



The model with Pq = 2 seconds (see §5 and Figs. 1, 2, 3, 
and 5) serves to illustrate the effects of viscous dissipation most 
clearly. Without including dissipation, at 290 ms after bounce 
there is a well-defined gain region with Rg ~ 115 km and a 
shock radius of ^ 275 km. The model is not exploding. It has 
timescale profiles very much like those in Fig. 4 for the model 
with Pq = 8 s with TAdv < Tq^ in the gain region (see §5.1). In 
this late post-bounce epoch we restart the calculation including 
viscous dissipation, ^'mri with a = 0.1 (see eq. 19). The reason 
for restarting the simulation with dissipation nearly 300 ms af- 
ter bounce instead of just 100 or 200 ms after bounce is that by 
this time in this progenitor the Si shell has fallen through the 
shock and the dynamical phase of shock expansion associated 
with this event is over. Because the system has settled into a 
steady-state accretion flow it is more easily diagnosed as we add 
heating by viscous dissipation. For this model, we include ^'mri 
only in the gain region above the PNS surface for r > 120 km. 
Fi gure 6 shows the relative contributions to <7tot — 4v + '7mri 
(solid line) at a time 340 ms after bounce. The heating rate 
from viscous dissipation (^mri), the total neutrino heating rate 
{Qv), and the contributions to from Ug and i>e neutrinos are 
shown. 

Complementary to Fig. 6, but more germane to the issue of 

explosion, is Fig. 7, in which we show all of the timescales 
presented in Fig. 4 and discussed in §5.1, but include r^^^, = 
{P / p) / qMRi (short dashed line). We distinguish between r^^ = 
iP/p)/q^ (dot-long dashed fine) and Tg^^ = (P/p)/(gV + ?MRi) 
(soUd line), where the former contains contributions only from 
neutrino heating and the latter includes both neutrino heating 
and viscous dissipation of shear energy. As shown clearly in 
Fig. 7, viscous dissipation is sufficient to make Tqj^j less than 
TAdv As per the discussion in §5.1 this is suggestive of a con- 
dition for a thermally-driven explosion. In this simulation and 
others, we find that if the condition TAdv/TV/TOT > 1 is met just 
outside Rg, the pressure increase is communicated to the rest 
of the gain region on a dynamical timescale and the shock ex- 
pands. That the number on the right-hand-side of this inequality 
is so close to unity is surprising considering the ambiguities in 
defining TAdv and Tg^Qj. Nevertheless, we find that it is accu- 
rate. The region in Fig. 7 where TAdv/^^ > 1 grows in time. 
The gain region expands as the shock moves outward. When 
the thermal energy of the gas exceeds the gravitational bind- 
ing energy (and the escape temperature obtains) an explosion 
develops (eq. 20). 

The explosion is not driven by viscous dissipation alone. It 
is triggered by the combined action of ^mri and q^. In fact, just 
^ 35 milliseconds after the condition rAdv/T^roi > 1 is niet in 
Fig. 7, TAdv has increased enough so that Tq^ < TAdv That is, 
neutrino heating alone is then sufficient to drive the explosion. 
In these ^ 35 milliseconds quRi deposits only ^ 2.1 x lO'*^ erg 
in the gain region. The neutrinos account for nearly five times 
that (Fig. 6). This illustrates explicitly that any heating source 
that contributes an additional ~ 20 - 50% to the net neutrino 
heating rate should be sufficient in most models to instigate ex- 
plosion (see §5.1 and Fig. 4). In the model described here it 
takes ~ 120 ms after the condition TAdv/^"^ > 1 is met for 
matter in the gain region to attain positive velocity. In this time 
period, 3.1 x 10"*^ erg of shear energy is dissipated, whereas 
~ 1.7 X 10^" erg is deposited by neutrinos. Because the explo- 
sion is initiated at late times (we did not add ^'mri until ^ 290 
ms after bounce), the extrapolated asymptotic energy of this ex- 
plosion is rather small. At the end of the calculation, 1 . 1 sec- 
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onds after bounce, taking rough account of the kinetic, thermal, 
and gravitational potential energies we estimate an asymptotic 
kinetic energy of 1.6 x 10^*^ erg. Including the energy gained 
by recombination to nuclei of the ~ 0.03 Mq with positive ve- 
locities at the end of the calculation adds ^ 5. 1 x 10^" erg, or a 
total asymptotic energy of ~ 6.7 x 10^*^ erg. Continued energy 
injection by neutrino heating will increase this estimate. The 
protoneutron star gravitational mass at the end of this calcula- 
tion is ~ 1.404 Mq. 

Including viscous dissipation only in the region r > 120 km is 
quite conservative. In fact, although there are strong stabilizing 
entropy gradients in the cooling region, because of the strong 
negative shear there, part of the cooling layer can be driven un- 
stable by the MRl. In addition, because this region does not 
encompass many pressure scale heights, we expect magneto- 
convective plumes to overshoot the gain region and enter the 
neutrino cooling layer. For these reasons we carried out two 
calculations analogous to the one just described, but which in- 
clude viscous heating in the cooling region (for r > 60 km and 
for r > 80 km). Both models explode at earlier times and with 
higher neutrino luminosities than the model with dissipation ex- 
terior to r > 120 km. Because the explosion is instigated ear- 
lier, these models leave behind neutron stars of smaller mass 
and yield explosions of higher energy. 

Similar numerical experiments may be carried out with dif- 
ferent initial rotation periods and as. For larger a and given 
Pq, the viscous time is shorter, shear energy is dissipated more 
rapidly (though roughly the same total energy is dissipated), the 
condition that TAdv/Tq^ > 1 is more easily obtained, and explo- 
sions take less time to develop. For smaller a, the threshold set 
by the condition TAdv/T^Tox > 1 "ot be reached or it may 
happen on a timescale too long to explore in the models pre- 
sented here.'^ For example, taking a = 0.01, qum/qv ~ 0.03 
in the gain region for the Pq = 2 s model just considered (see 
Fig. 6), TAdv/T'.jTOT < 1 throughout the gain region, and no ex- 
plosion is obtained. 

For shorter Pq (faster rotation), more shear energy is avail- 
able in the gain region and at a given a an explosion develops 
more quickly. In fact, the viscous heating rate at a given radius 
scales with fP and, all else being equal, if Pq is decreased from 
2 s to 4 s, the heating rate in the gain region decreases by a 
factor of ~ 6. Thus, for a given a, beyond a threshold Pq the 
difference between tajv and r^^ cannot be compensated for by 
including ^mri and no explosion obtains. However, this thresh- 
old depends on the extent of the region of viscous heating. 

It is perhaps artificial to begin including viscous dissipation 
so long after bounce. The purpose of the previous experiment 
was merely to clearly isolate the physics associated with shear 
energy dissipation. In reality, the magnetic field grows from an 
arbitrarily small seed field exponentially on a timescale ^ . 
We expect a saturation magnetic field to obtain and fully de- 
veloped magneto-convection to appear just 10-100 ms after 
bounce in the region exterior to the PNS. For this reason, we 
have run a number of models including ^mri after bounce, fol- 
lowing the full evolution. 

To this end, we conduct two tests with the model E15A (see 
§3.2, Table 1, and Figs. 1 and 2). In one case, we include vis- 
cous dissipation at all radii. In the second calculation we turn on 
viscous dissipation only when the generalized Solberg-H0iland 

Due to numerical limitations, including shock resolution and Courant-limited tii 



criterion indicates that the region is unstable (eq. 12). This latter 
case is probably more realistic within the context of the spheri- 
cal models presented in this paper. Both models employ a = 0. 1 
and do not transport angular momentum until after bounce. In 
both cases, we obtain explosions that develop on the relatively 
short timescale of ^ 50- 100 ms after bounce. 

The explosions in this model are qualitatively different from 
the explosion obtained by waiting until a well-defined gain re- 
gion forms and then adding viscous dissipation (Figs. 6 and 7). 
In non-exploding models, the bounce shock stalls immediately 
after formation, but then moves slowly outward in radius as the 
cooUng and gain regions are formed until the shock reaches 
/?sh ~ 200 km (about ^ 100-200 ms after bounce). Its outward 
progress is halted primarily by neutrino losses in the forming 
cooling region. Usually, the shock radius recedes slowly after 
reaching its maximum radius and sits for many hundreds of mil- 
liseconds at r ~ 150 km. Indeed, this evolutionary sequence is 
borne out in model E15A without viscous dissipation; no explo- 
sion occurs. Including viscous heating in E15A, the evolution is 
very different. Neutrino losses behind the shock (as R^h moves 
from ~80 km to ^200 km) are partially compensated by vis- 
cous heating. This extra energy deposition allows the shock to 
continue its outward progress, putting matter higher in the PNS 
gravitational well, and explosions are obtained. Because the 
shear in this model is rather large, the region behind the shock 
is unstable to the MRI and viscous dissipation operates. Stabi- 
lizing entropy gradients are insufficient to make N^+dfP'/dlnr 
positive (eq. 12). Because both the calculation with dissipation 
everywhere and the calculation respecting eq. (12) have vis- 
cous heating that partially compensates neutrino cooling as the 
shock is moving outward, the same qualitative evolution is ob- 
tained in both models. One important difference between them 
is that, on average, the model with dissipation throughout the 
profile has more heating behind the shock than the model that 
employs the generalized Solberg-H0land criterion. This means 
that the explosion develops faster and leaves behind a less mas- 
sive neutron star. Although we are able to follow only the first 
~ 200 ms after explosion (^ 500 ms after bounce), the differ- 
ence between the PNS mass at the end of our simulations is 
roughly ^ 0. 15 Mq, the larger baryonic mass being ^ 1 .4 Mq. 

Figure 8 shows the angular momentum evolution for the 
EISA model with dissipation at all radii as a function of the en- 
closed mass. For reference, the dotted line shows jKep = r^^Kep 
at bounce (t = 0). The profile j(M) is shown at bounce and 22, 
100, 200, and 350 ms after bounce. The transport of angular 
momentum is clear. For high viscosity, f2 can decrease appre- 
ciably throughout the PNS core during explosion, providing a 
mechanism for at least partial spindown of the young PNS. Fig- 
ure 8 can be compared with the three-dimensional simulations 
of Fryer & Warren (2004), who found significant angular mo- 
mentum transport in the deep core of their rotating progenitors 
(see their Figs. 9 and 10). At M ~ 0.5 Mq they find that ; 
decreases by a factor of ^ 4 in ^ 200 ms. We find a similar 
magnitude change in j at this M, implying that the magnitude 
of the intrinsic angular momentum transport in the simulations 
of Fryer & Warren (2004) is comparable to the transport caused 
here by our expUcit viscosity. 

Similar to the E15A models just described, we computed 

the evolution of Po = 1.25, 2, 3, 4, and 6 s models with the 
11 Mq progenitor including viscous dissipation with a = 0.1. 

stepping, it is difficult to push our calculations much beyond 1 s after bounce. 
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Like the E15A models, we calculate a set of models with dis- 
sipation everywhere at all radii and a set of models respecting 
the generalized Solberg-H0iland stability criterion. All models 
with Po < 4 s explode energetically before a well-defined gain 
region is established. The model with Pq = 4 and dissipation 
everywhere also explodes, but on a longer timescale (~ 200 
ms after bounce). The model with Pq = 6 s does not explode. 
The Pq = 4 s model with dissipation only in regions where the 
MRI is unstable also does not explode. In this model, dur- 
ing the post-bounce epoch, quRi was not large enough to force 
< TAdv in the gain region. In addition, just after bounce 
the shear is not strong enough to drive a large portion of the 
neutrino cooling region unstable to the MRI. For this reason, 
explosions qualitatively similar to the E15A models do not ob- 
tain. In summary, our calculations and estimates of the energy 
deposition rate of viscous dissipation show that for initial spin 
periods in the range Po < 4-6 s in our 1 1 Mq progenitor, and 
employing the angular momentum profile dictated by equation 
(4), energetic supernova explosions can be obtained for a ~ 0. 1 
with dissipation everywhere. We therefore quote an initial spin 
period required to significantly modify the dynamics and out- 
come of one-dimensional simulations of core-collapse super- 
novae of Po ^ 5 s. More detailed predictions must await multi- 
dimensional MHD simulations. 

Typical results for the explosions we obtain are shown in 
Figs. 9 and 10. The former shows snapshots of the velocity 
profile at various times during the development of explosion in 
our model with Pq = 2 s. This particular model included viscous 
dissipation at all radii, but the results for this initial spin period 
are very much like those with dissipation only in regions where 
the generalized Solberg-H0iland stability criterion indicates in- 
stabiUty. The shock stalls and an accretion phase begins, but 
this model never reaches a steady state with a well-defined gain 
region. The constant injection of thermal energy by viscous 
dissipation keeps the pressure behind the shock relatively high 
and it moves out slowly, putting post-shock material higher in 
the PNS's potential well. An explosion develops on a timescale 
of just ~ 100 ms after bounce. In Figure 10 we show the full 
history of the mass element trajectories in radius as a function 
of time in the model with Pq = 4 s and dissipation everywhere. 
This explosion takes considerably more time to develop. In this 
model it is primarily the viscous heating in the cooling region 
that instigates explosion after a well-defined gain region has de- 
veloped. The heavy solid lines show the trajectories of the 1 .0, 
1.2, 1.3, and 1.4 Mq mass elements. The simulation is very 
well resolved; only every tenth mass zone is shown. The neu- 
tron star mass at the end of the calculation is ~ 1 .35 Mq. 

As an aside, in Fig. 1 1 we provide two snapshots of the ve- 
locity profile (solid lines) in a representative exploding model. 
For comparison, we also include the sound speed profile. The 
two snapshots are separated by ~ 100 ms in time and the shock 
moves from approximately 750 km to 2000 km. The large dot 
at '~850 km marks the sonic point and the emergence of a tran- 
sonic PNS wind (Duncan et al. 1986; Woosley et al. 1994; 
Takahashi, Witti, & Janka 1994; Burrows, Hayes, & Fryxell 
1995). Profiles of entropy, mass loss rate, and electron fraction 
all show that a quasi- steady- state solution, quantitatively simi- 
lar to that expected from the models of Thompson et al. (2001) 
and the analytic estimates of Qian & Woosley (1996) is ob- 
tained as the PNS begins the Kelvin-Helmholtz cooling epoch. 
In all of our successful exploding models, a transonic wind is 
formed. At the end of this calculation (Pq = 2 s, a = 0.1), 100 



ms after the second snapshot in Fig. 11, the mechanical power 
of the wind is > 10^" erg s"' . If strong transonic PNS winds are 
a generic feature of the core-collapse supernova phenomenon, 
as implied by our calculations, this conclusion may have im- 
plications for late-time fallback and the mass cut (Woosley & 
Weaver 1995). 

8. SUMMARY, CONCLUSIONS, & IMPLICATIONS 

In this paper we forward the idea that some of the gravita- 
tional binding energy of core collapse may be trapped in differ- 
ential rotation and that this free energy can be tapped on a vis- 
cous timescale (^ 1 s in the gain region) and deposited as ther- 
mal energy in the gas. We identify magnetic stresses generated 
by the magnetorotational instability as the most likely origin of 
this viscous energy dissipation. We have constructed models of 
core collapse, adding rotation and viscosity to our existing al- 
gorithm for spectral, multi-angle radiation hydrodynamics. We 
find that viscous dissipation can enhance the heating in the gain 
region sufficiently to yield explosions in models that would oth- 
erwise fail. Initial spin periods of the roughly solid-body iron 
core should be Pq < 5 seconds for this mechanism to be robust. 

We emphasize several important points about this result. 
First, our simulations are one-dimensional and as such do not 
include the convective motions that have been shown to at least 
aid, if not enable, explosions (Herant et al. 1994; Burrows, 
Hayes, & Fryxell 1995; Janka & MuUer 1996; Fryer et al. 1999; 
Janka et al. 2002; Buras et al. 2003). We expect that multi- 
dimensional models with slower initial progenitor rotation pe- 
riods, including the effects of viscous dissipation consistently, 
will also yield explosions. Second, the effects of dissipation 
depend on where, when, and how quickly dissipation acts. Im- 
proved understanding of these issues will require MHD simu- 
lations and further progenitor modeling. Third, the values of 
Pq we require to significantly effect the dynamics of core col- 
lapse are fully consistent with the rotating progenitor models of 
Heger et al. (2000). However, there is considerable uncertainly 
in these progenitor rotation rates. In particular, the models of 
Heger et al. (2003) and Heger et al. (2004), which attempt to 
include angular momentum transport by magnetic torques con- 
sistently in one dimension, find much slower rotation rates for 
the iron core just prior to collapse. For low mass progenitors in 
the range 1 2 - 1 5 M© , the rotation periods obtained by Heger et 
al. (2004) are larger than our Pq < 5 s fiducial bound. Heger et 
al. (2004) do, however, find that the rotation period decreases 
with increasing progenitor mass. Thus, without a detailed study 
beyond the scope of this paper, we cannot make a definite state- 
ment as to whether rotational effects are important in their core- 
coUapse progenitor models. 

Finally, our limit on Pq is, at face value, inconsistent with the 
observed distribution of spin periods of young pulsars, which 
are in the range of many tens to many hundreds of millisec- 
onds (Kaspi & Helfand 2002). There is, of course, a danger in 
such comparisons since it is very difficult to establish the true 
birth periods of neutron stars observationally. It is also unclear 
if the spindown of very young neutron stars is well-understood 
theoretically as there are a number of ways to remove angular 
momentum from the neutron star during the explosion phase 
(as in the models presented here with viscous dissipation) or 
during the Kelvin-Helmholtz cooling epoch (§2). Given these 
uncertainties, we reiterate that the range of rotation rates con- 
sidered here are expected theoretically in models of magnetar 
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formation (Duncan & Thompson 1992; Thompson & Duncan 
1993; Kouveliotou et al. 1999), a possibly large subset of all 
supernovae. 

Our work has highlighted two primary impacts of rotation 
on the mechanism of core-collapse supernovae: 1) Most impor- 
tantly, rotational kinetic energy tapped by viscous processes can 
help drive explosions (§7), and 2) Rotation lowers the effective 
gravity in the core, increasing the radius of the stalled shock 
and the size of the gain region (§5, Fig. 3). Since ejection is 
inhibited by the deep potential well, this helps facilitate explo- 
sion. In addition to these, there are several multi-dimensional 
effects of rotation on core-collapse that remain to be fully ex- 
plored (Burrows, Ott, & Meakin 2003): 1) Rotation might gen- 
erate vortices that can dredge up heat from below the neutri- 
nospheres and thereby enhance the driving neutrino luminosi- 
ties. 2) Rotation can result in a large pole-to-equator anisotropy 
in the mass accretion rate after bounce due to the centrifugal 
barrier along the poles. Finally, 3) rotation results in pole-to- 
equator anisotropics in the driving neutrino flux that have never 
before been accurately calculated. Since the polar neutrino flux 
should be enhanced and the polar mass flux should be lowered 
(point 2), explosion along the poles with bipolar morphology 
seems that much more likely (Kotake, Yamada, & Sato 2003; 
Burrows & Goshy 1993; Burrows, Ott, & Meakin 2003). Thus, 
such bipolarity, and the consequent optical polarization of the 
debris, are not exclusive signatures of MHD-driven explosions 
(Akiyama et al. 2003; Symbalisty 1984), but may be a natural 
consequence of any rapidly rotating models. 

In addition to the effects of rotation on the supernova mech- 
anism, we have also shown how the observational signature 
of core collapse might be affected by rotation by identify- 
ing the systematics in the neutrino luminosity and average en- 
ergy as a function of time, initial rotation rate, and progenitor 
(§5, Fig. 2). If the neutrino radiation from the next galactic 
supernova is observed in modern neutrino detectors like Su- 
perKamiokande and SNO, we may hope to distinguish rapidly 
rotating core collapse via its neutrino signature. In addition to 
the difference highlighted in Fig. 2 we note that rapidly rotat- 
ing models may leave an appreciable amount of shear energy 
(^ 10^^ erg) trapped in the inner core, which is stable to the 
MRI and convection. The timescale to dissipate this energy 



is uncertain, but could be ~ 1-10 seconds or longer (eq. 21; 
Spruit 2002) as the PNS cools. The timescale for dissipa- 
tion of this energy will likely differ from the Kelvin-Helmholtz 
timescale (Burrows & Lattimer 1986; Pons et al. 1999). Thus, 
rapidly rotating PNSs may have a unique long-timescale neu- 
trino signature. 

Finally, we have assessed the relative role of neutrino vis- 
cosity, the turbulent viscosity caused by the magnetorotational 
instabiUty, and the viscosity of purely hydrodynamic convec- 
tion. We find that, for the rotation rates considered here, the 
MRI and convection may be comparable in importance, but that 
both dominate microscopic neutrino viscosity everywhere and 
at aU times during the collapse and accretion epochs. Turbulent 
magneto-convection should be generic for all modestly rotating 
progenitors and of critical importance to the dynamics of mod- 
els with fairly rapid rotation, Po < 5 s (Duncan & Thompson 
1992; Thompson & Duncan 1993). Thus, in models with these 
rotation rates, magnetohydrodynamic effects should not be ig- 
nored. The simplicity of the algorithm employed here, which 
assumes spherical symmetry in a system that so clearly violates 
it and an a viscosity, makes it clear that more work must be 
done to ascertain the full effects of MHD in rotating core col- 
lapse. 
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Table 1 

Rotating Core-Collapse Without Dissipation 



Name 


^Rot 


/3r 


oRot 

Psoo 




El^^ (1051 erg) 


^0= 1.25 s 


0.0035 


0.059 


0.070 


21.8 


32.9 


Po = 2.0 s 


0.0014 


0.028 


0.039 


11.4 


25.0 


Po = 3.0s 


6.0 X 10-4 


0.012 


0.020 


4.98 


13.9 


Po = 4.0s 


3.4 X 10-4 


0.0065 


0.012 


2.46 


8.39 


/^o = 5.0s 


2.2 X 10-4 


0.0042 


0.0077 


1.61 


5.56 


Po = 8.0 s 


8.5 X 10-5 


0.0017 


0.0031 


0.638 


2.26 


E15A 


0.0020 


0.053 


0.067 


20.6 


49.7 



0^°'^ is the ratio of the total rotational energy to the magnitude of the total gravitational energy. All "Pq" models have Rn = 1000 km (see eq. 4). 
Subscripts /, b, and 500 imply initial, at bounce, and 500 milliseconds after bounce, respectively. 
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Fig. 1 . — Initial angular velocity (Q, rad s"') as a function of radius for the models considered in this paper The soHd lines show the Q(r) profile imposed on the 
1 1 Mq model of Woosley & Weaver (1995) with = lOOO km and Pq = 1-25, 2, 3, 4, 5, and 8 seconds (see eq. 4). The dashed line shows the angular velocity 
profile for the model E15A from Heger et al. (2000). For simplicity, we refer to the Woosley & Weaver (1995) models by their initial "Pq" and the Heger et al. (2000) 
model as "E15A". 

Fig. 2. — Neutrino luminosity (L,y, left panels) in 10^^ erg s"' and average neutrino energy ((e^), right panels) in MeV as a function of time after bounce for the 
rotating 11 Mq models with Fu = 1.25, 2, 3, 4, 5, and 8 s (solid lines). In all cases, at 200 ms after bounce, the lowest Li, and (s^) correspond to the fastest initial 
rotation. Note that this hierarchy is preserved at all times and in all panels except in the case of L^^ 400 ms after bounce and at breakout. Although they cannot be 
directly compared with the 11 Mq models, the results for the model E15A from Heger et al. (2000) are also shown (dashed Unes, see Fig. 1) The Ue breakout burst 
at t = (not shown here) is largest for the fastest rotator, with L^f^ ~ 3. 1 x lO'^ erg s"' . For both E15A and the model with Pq = 2 s, L^f" ~ 2.6 x 10" erg s"' . 
For the slowest rotator we recover the non-rotating case and L^f^ ~ 2.4 x lO'' erg s~' . See §5 for details. 

Fig. 3. — Mass accretion rate (M [Mq s~'], upper left panel), shear (dCl/dlnr) and CI (10^ rad s~\ upper right panel), putative saturation magnetic field strength 
(logjQ[Bsat G], lower left panel), and viscous dissipation rate (^mri [erg g~^ s"'], lower right panel) versus radius in km for rotating models with Pq = 1.25 (dotted), 
2 (dashed), 3 (dot-dashed), and 8 s (soUd) at r = 105 ms after bounce. These models did not include viscous dissipation. We show ^mri (with a = 0. 1 in ^mri as in 
eq. 15) only for comparison with typical neutrino heating rates in models of supemovae (q^, soUd box). For clarity of presentation a different radial scale is used 
for each panel. 

Fig. 4. — The neutrino heating timescale (th,^ = (Pj p)/Hi,, dot-dashed line), neutrino cooling timescale (tq^ = (P/ p)ICi,, dotted line), net neutrino heating 
timescale (r^^ - {P/ p)/qu, solid line), and the advection timescale (TAdv = H/Vr, dashed line) in seconds, as a function of radius r (km) at a time 130ms after 
bounce. Rv^ and ffg mark the ve neutrinosphere (for sv^ ~ 12 MeV) and the gain radius, respectively. 

Fig. 5. — Left Panel - Microscopic shear viscosity of electron neutrinos (solid line), anti-electron neutrinos (long dashed line), and muon neutrinos (short dashed 
line) in the deep interior of a PNS ~105ms after bounce. For comparison, we show the effective viscosity due to the MRI (taking a = 0.1; dotted Une), which 
dominates the total neutrino viscosity. Right Panel - Damping rate due to neutrino interactions (F^) and the optimal growth rate of the MRI (Fmri ~ H) in the 
semi-transparent regime. The Brunt- Vaisala frequency is also shown. 

Fig. 6. — Energy deposition rate as a function of radius in units of 100 MeV baryon"' s"' for the rotating model with Po = 2 seconds, including viscous dissipation 
via the MRI with a = 0.1. Total energy deposition rate (^tot = 4v+4msi, soUd Une), viscous energy deposition rate (^mri, eq- 19, dot-dashed line), and neutrino 
deposition rate qi, =Hu-Cu, long dashed Une) are shown. For comparison, individual contributions to qi, from Ve (short dashed Une) and Pe (dotted Une) are also 
included. 

Fig. 7. — The neutrino heating timescale (th^ = (P/ p)lHi^, dot-short dashed line), neutrino cooling timescale (tq^ = (P j p)jCv, dotted line), net neutrino 
heating timescale {Tq-^ — {P/ p)/q,j. dot-long dashed line), the advection timescale (taiJv = H /Vr, long dashed line), the viscous heating timescale for the MRI 
("TvlRl = (P / p) I qyiKl, short dashed line), and the total heating timescale including viscous dissipation (t^jqi ~ I pMiiv +^mri), solid line) in seconds, as a function 
of radius r (km) at a time 340 ms after bounce in our 1 1 Mq progenitor with Pq = 2 s. Ru,, and Pg mark the neutrinosphere and gain radius, respectively. The 
shock radius is ~ 275 km. Compare with Figs. 4 and 6. Note that exterior to Pg, as a result of ^mri, TAdv is greater than (see §5.1 and §7). 

Fig. 8. — Angular momentum (j, lO'^ cm- s"' ) as a function of enclosed mass in the model E15A including viscous dissipation with a = 0.1 at bounce (solid line, 
t = 0), 22 ms (long dashed line), 100 (dot-short dashed line), 200 (short dashed line), and 350 ms (dot-long dashed line) after bounce. For reference the Keplerian 
angular momentum profile 7'Kep = '^f^Kep at boimce (t = 0, dotted Une) is also shown. 

Fig. 9. — Radial velocity as a function of radius for various snapshots in time during the development of explosion in the 11 Mq progenitor with Pq = 2 s, 
including viscous dissipation after bounce with a = 0.1. 

Fig. 10. — Radial position of mass elements as a function of time in a model with Pq = 4 s and a = 0. 1. Viscous dissipation was assumed to occur throughout the 
radial profile. The heavy solid lines show the trajectories of the 1.0, 1.2, 1.3, and 1.4 Mq mass elements. At the end of the simulation, the mass cut is at ~ 1.35 
Mq . Note that for clarity only every tenth mass zone is shown in this figure. 

Fig. 1 1 . — Matter velocity and adiabatic sound speed at two snapshots during explosion in the model with Pq = 2 s (compare with Fig. 9). Note the development of 
the sonic point on a timescale of ~ 100 ms, indicating the presence of a transonic wind. This neutrino-driven super-sonic outflow is generic to aU of our exploding 
models. 
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